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Correlation induced switching of local spatial charge distribution in two-level system 



O 

(N 

C 

0\ 



i 

c 

o 
o 



> 

(N 

o 



X 



P. I. ArseyevQ N. S. MaslovaQ and V. N. MantseviclH 
P.N. Lebedev Physical institute of RAS, 
119991, Moscow, Russia 

Moscow State University, 
Department of Physics, 119991 Moscow, Russia 
(Dated: January 28, 2013) 

We present theoretical investigation of spatial charge distribution in the two-level system with 
strong Coulomb correlations by means of Heisenberg equations analysis for localized states total elec- 
tron filling numbers taking into account pair correlations of local electron density. It was found that 
tunneling current through nanometer scale structure with strongly coupled localized states causes 
Coulomb correlations induced spatial redistribution of localized charges. Conditions for inverse oc- 
cupation of two-level system in particular range of applied bias caused by Coulomb correlations have 
been revealed. We also discuss possibility of charge manipulation in the proposed system. 
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I. INTRODUCTION 

Investigation of tunneling properties of interacting im- 
purity complexes in the presence of Coulomb correlations 
is one of the most important problems in the physics 
of nanostructures. Tunneling current changes localized 
states electron filling numbers as a result-the spectrum 
and electron density of states are also modified due to 
Coulomb interaction of localized electrons. Moreover the 
charge distribution in the vicinity of such complexes can 
be tuned by changing the parameters of the tunneling 
contact. Self-consistent approach based on Keldysh di- 
agram technique have been successfully used to analyze 
non-equilibrium effects and tunneling current spectra in 
the system of two weakly coupled impurities (when cou- 
pling between impurities is smaller than tunneling rates 
between energy levels and tunneling contact leads) in the 
presence of Coulomb interaction [ij. In the mean-field 
approximation for mixed valence regime the dependence 
of electron filling numbers on applied bias voltage and 
the behaviour of tunneling current spectra have been an- 
alyzed in 0- 

Electron transport even through a single impurity in 
the Coulomb blockade and the Kondo regime [H have 
been studied experimentally and is up till now under 
theoretical investigation As tunneling coupling 

is not negligable the impurity charge is not the discrete 
value and one has to deal with impurity electron fill- 
ing numbers (which now are continuous variables) de- 
termined from kinetic equations. 



Analyzing non-equilibrium tunneling processes 
through coupled impurities one can reveal switching 
on and off of magnetic regime (electron filling numbers 
in the localized states for opposite spins are equal) on 
each impurity atom at particular range of applied bias 
voltage [2|- 

In the present work we consider the opposite case when 
coupling between localized electron states strongly ex- 
ceeds tunneling transfer rates. This situation can be ex- 
perimentally realized when several impurity atoms or sur- 
face defects are situated at the neighboring lattice sites, 
so coupling between their electronic states can strongly 
exceeds the interaction of this localized states with con- 
tinuous spectrum Figjl] Another possible realization 
is two interacting quantum dots on the sample surface 
weakly connected with the bulk states. Such systems 
can be described by the model including several electron 
levels with Coulomb interaction between localized elec- 
trons. If the distance between impurities is smaller than 
localization radius, strong enough correlation effects arise 
which modify the spectrum of the whole complex. Elec- 
tronic structure of such complexes can be tuned as by 
external electric field which changes the values of sin- 
gle particle levels as by electron correlations of localized 
electronic states. One can expect that tunneling current 
induces non-equilibrium spatial redistribution of local- 
ized charges and gives possibility of local charge den- 
sity manipulation strongly influenced by Coulomb cor- 
relations. In some sense these effects are similar to the 
"co-tunneling" observed in [TTJ] , 12 1. Moreover Coulomb 
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interaction of localized electrons can be responsible for 
inverse occupation of localized electron states. These ef- 
fects can be clearly seen when single electron levels have 
different spatial symmetry. 

To understand such correlation induced "charge" 
switching it's sufficient to analyze Heisenberg equations 
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for localized states total electron filling numbers tak- 
ing into account pair correlations of local electron den- 
sity [lH. If one is interested in kinetic properties and 
changes of local charge density for the applied bias range 
higher than the value of energy levels tunneling broad- 
ening modification of initial density of states due to the 
Kondo effect can be neglected. In this case for the finite 
number of localized electron levels one can obtain closed 
system of equations for electron filling numbers and their 
higher order correlations. 



with elctron levels Operators c^ p ^/c fe ( p ) correspond 
to electrons creation/annihilation in the continuous spec- 
trum states k(p). rn a = c^c^-two-level system electron 
filling numbers, where operator Ci„ destroys electron with 

spin a on the energy level U^ a is the on-site Coulomb 
repulsion of localized electrons. 

Tunneling current through the two-level system can 
be written in the terms of electron creation/annihilation 
operators as: 



II. THE SUGGESTED MODEL 



We shall analyze tunneling through the two-level sys- 
tem with Coulomb interaction Fig[T] The model system 
can be described by the Hamiltonian H . 

a 




Sample 




I = ha = ^ n k(7 = t k (< c+ a c ia > - < cf a c ka >)(2) 



Let us consider h — 1 elsewhere, so motion equation 
for the electron operators product ct g .Ci <r can be written 



; ® C ka Cia 
' dt 



+ Uij(nj a + rij-a) ■ c i a CicT - tk ■ (n irT - f k ) + 

k' #fc i^j 



where 



(4) 



FIG. 1: a). Energy diagram of two- level system and 
b). Schematic spatial diagram of experimental realization. 
Coulomb energy U12 correspond to the interaction between 
electrons on different energy levels. 



ia ka per 



hi a 



(1) 



Now let us also consider that n| CT = n itT . 

Neglecting changes of electron spectrum and local den- 
sity of states in the tunneling contact leads due to the 
tunneling current flowing we shall uncouple conduction 
and two-level system electron filling numbers. After sum- 
mation over k one can get an equation which describe 
tunneling current in the presented two-level system: 



I k ia + Ik 



I,-2n 



(5) 



Indices k and p label continuous spectrum states in 
the left (sample) and right (tip) leads of tunneling con- 
tact respectively. t k ( p y tunneling transfer amplitudes be- 
tween continuous spectrum states and two-level system 



Where expression for tunneling current I k 2 a can be 
obtained by changing indexes I o- 2 in equation for tun- 
neling current I k x a which has the form: 
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hia = T fc • {(n la ) + y^(ct.c icr ) - 

- {(1 - ni_ CT )(l - n 2 _ CT )(l - n 2(T )} ■ /fe(ei) - 

- (ni_ CT (l - n 2 _ CT )(l - n 2cr )} • / fc (ei + E/n) - 

- E< n 2 CT ' (1 - n 2 _„'){l - ■ / fc (ei + C/12) - 

- 5Z^ 1 - ct71 2 CT '( 1 -«2-<7')> ' + Uu+U 12 ) ~ 
(j 

- {n2an 2 - a {l - Tll-a)) ■ fk{si + 2J7i 2 ) - 

- (ni- a n2- a n2 a ) ■ fk(ei + Uu + 2f7i 2 )} + 

(1 - -n 2 -a)(l -«2<r)\ . 



{( 



£1 — £fc 



- n 2 - CT )(l - n 2<T ) , 
ei + C/11 - e fe 

+ , Sg-"v( 1 -"l^)( 1 - Tt 2V) i | 

£l + t/l2-£fc 



+ 



' e 1 + Un + Ui2-s k 
,712-^712^(1 - ni- a ) 



£1 + 2i7i 2 - £fc 



) + ( 



"l-g-»2-g-»-2cr 
El + + 2t/i2 - £fe 



)} 



(0) 



We shall further neglect terms tkcf^Cja and tkC~l a c k > a 
in expression [6] as they correspond to the next order per- 
turbation theory by the parameter ^g- . Relaxation rates 

Tfc(p) = 7r • • i^o are determined by electron tunneling 
transitions from two-level system to the leads k (sample) 
and p (tip) continuum states, fo-continuous spectrum 
density of states. Equations for filling numbers ni a 77,20- 
can be found from the conditions: 



dn lcr 

at 

dn 2t y 



Ikln + Iplcr — 



Ik2o 



lp2<7 







(7) 



where tunneling current I pa can be easily determined 
from Iko- by changing indexes fc «-» p 

We shall analyze the situation when Coulomb energy 
values are large and condition E/y >> £j/j can be taken 
into account. It means that if one have to calculate tun- 
neling current through such system it is necessary to find 
all pair filling numbers correlators in the energy range 
Si + Uij. So we retain the terms containing fk( p ) (£« + Uj) 
and neglect all high orders correlators and pair correla- 
tors which contain / fe ( p ) (si + Uij + Uki). We consider the 
paramagnetic situation 71^ = rii^ a . 

Pair filling numbers correlators can be found in the 
following way: 



dn ia n ja < dn ia dn j(T > 

<— 5T-> = <-5r^> + (-^-"*r> (8) 

Let us introduce tunneling filling numbers n T (si), 
T (si + Uij) and nf, which have the form: 



n T (£i) = 



n T (£i + Uij) = 



where 



Tkfk(ei)+T p f p (ei) 

Tkfk(e t + U lJ ) + T p f p (e l + U lJ ) 
T k +T p 



(9) 



fkij — fk( £ i) — fk( e i + Uij) 



(10) 



As we consider that m a — ni^ ai let us also consider 
(uiaTija) = (nio-rij-o-). So a system of equations for 
pair correlators Ku = (nxo-nx-^), K22 = {n2<y n>2-a) an d 
K12 = {ni a ri2a) for large Coulomb energies Uij » s^u 
has the form: 



an 


a.12 


ai 3 \ 
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a-22 


(123 




K- 
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a-32 


033/ 







= F 



(11) 



where 



an 
ai3 
ai2 
a 2 2 



023 = 1 
a 2 i = 
2tt t (£ 1 + U u ) 
2n T (e 2 + U22) 



(12) 



031 

032 
033 



^« T (£2 + U12) 



= 1 



1 T 1 

-n ( £l 



U12 



^n T (s 2 + Ui2) 



-n T (ei+U 12 ) 



(13) 



and 



F 



n T (ei + Uu) ■ ni a 
n T (e 2 + U22) ■ n 2a - 
2 ' 



)n T {ei + U12) ■ n 2<y + tU t (e 2 + U 12 ) ■ 77i CT 



(14) 



Equations which determine two-level system filling 
numbers immediately follows from the system 1111 
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bias voltage, eV bias voltage, eV bias voltage, eV 



FIG. 2: Two-level system filling numbers a).-c). and tunneling current d).-f). as a function of applied bias voltage in the case 
when both energy levels are situated above the sample Fermi level. Parameters ei = 0.6, £2 = 0.3, U12 = 1.0, Un = 1.4, 
U22 = 1.5 are the same for all the figures. a),d).r fc = 0.01, F p = 0.01; b),e)T k = 0.05, T p = 0.01; c),f).F fe = 0.01, T p = 0.03. 



n 2o 



(1 



n n 
+ 2-K 12 

= n T { £l ) 

~T 
n 22 



) + n 2o 



2nf a - K22 



(n 1 ( £l ) - 2n J (ei + U 12 )) + 



ci = £2 = £ and consequently Uij = U. In this case 
filling numbers have the form: 



(-n T ( £l ) + n T ( £l + U n ) + n T (£i + U 12 j) = 



(1 
2-K 12 

n T {e 2 ) 



) + ni 



2ni; 



^ii-(n T (e 2 )-2n T (e 2 + C/ 12 )) + 



(18) 



(-n T {e 2 ) + n T (e 2 + U 22 ) + n T (e 2 + U 12 j) = 



And finally expression for tunneling current has the 
form: 



hia = T fc • {(n la ) - (1 - (ni„) 
■ fk{ei)-{(n la )-2K 12 ) 
- 2 • ((n 2a ) - K 12 - K 22 ) 



- 2(n 2rT ) + K 22 + 2K 12 ) ■ 
fk(ei + U 12 ) (16) 



Let us also mention two extreme cases. The first one 
when all Coulomb energies are extremely large U%j — > 00. 
In this situation expressions for filling numbers will have 
the following form: 



>i)'(l 



n 2o 



(l+n^(e 1 ))-(l + n^(e 2 )) 
n T {e 2 ) ■ (1 - 7 



4 • n T (si) 



>2) 



(1 + n T { £l )) ■ (1 + n T {e 2 )) - 4 • n T (ei) • n T (e 2 ) 



(17) 

And the second one is when energy levels are gener- 
ated, for example due to the orbital quantum number 



(15) 
III. 



MAIN RESULTS AND DISCUSSION 



The behaviour of non-equilibrium electron filling num- 
bers with changing of applied bias and tunneling con- 
ductivity characteristics obtained from equations ((5]) and 
([TT ]> -(fTB j) are depicted on FigHFigH] 

We consider different experimental realizations: both 
energy levels are situated above the sample Fermi level 
(Figj2]); both levels below sample Fermi level (Fig(4j) and 
one of the energy levels is located above the Fermi level 
and another one below the Fermi level (Figj3j) . From the 
obtained results one can clearly see charge redistribution 
between two electron states with changing of applied bias 
voltage (Fig BIJ) . 

When both levels are situated above (FigJ^J) or below 
(Fig|4j the sample Fermi level one can clearly reveal two 
possibilities for charge distribution in the two-level sys- 
tem. The first one corresponds to the case when local 
charge is mostly accumulated on the lower electron level 
Tlx < n 2 (e 2 < eV < E\, e 2 + U\ 2 < eV < e% + U\ 2 and 
e 2 + U 22 < eV < £1 + U u on FigH and Figg}. The 
second one deals with the case when charge is localized 
on both levels equally n\ = n 2 [e\ < eV < £ 2 + U\ 2 , 
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FIG. 3: Two-level system filling numbers a).-c). and tunneling current d).-f). as a function of applied bias voltage in the 
case when one energy level is situated above and another one below the sample Fermi level. Parameters ei = 0.2, £2 = —0.3, 
U12 = 1.0, Uu = 1.4, U 22 = 1.7 are the same for all the figures. a),d).r fc = 0.01, F p = 0.01; b),e)T* = 0.05, T p = 0.01; 

c),f).r* = 0.01, r p = 0.03. 



£1 + U12 < eV < e 2 + U22 and e 1 + Uu < eV on FigfJ 
and FigHJ. 

Coulomb correlation induced sudden jumps down and 
up of each level electron filling numbers at certain values 
of applied bias are clearly seen. 

So if electron states have essentially different symmetry 
one can expect charge accumulation in various spatial ar- 
eas and thus the possibility of local charge manipulation 
appears. 

When both electron energies are situated below the 
sample Fermi level upper electron state become empty 
(ni = 0) for two ranges of applied bias voltage (£2 < 
eV < £1 and e 2 + U 12 < eV < £1 + Ui 2 ) (FigHJ) - 

Described peculiarities take place for all the ratios be- 
tween tunneling transfer rates Tk and T p . 

The other interesting effect is the possibility of inverse 
occupation of the two-level system due to Coulomb in- 
teraction in special range of applied bias (Figj3]). In the 
absence of Coulomb interaction difference of electron fill- 
ing numbers is determined by electron tunneling rates 
m — n,2 ~ 7/ci7p2 — 7pi7fc2- So without Coulomb inter- 
action, for 7fc(p)i = 7fe( P )2i difference of the two levels 
occupation numbers turns to zero. Coulomb interaction 
of localized electrons in the two-level system results in 
inverse occupation of two levels at the high range of ap- 
plied bias voltage. This situation is clearly demonstrated 
on the FigH 

It is clearly evident (Fig{3]) that when applied bias 
doesn't exceed value £1 + U12 all the charge is localized 
on the lower energy level (m = 0). With the increas- 
ing of applied bias inverse occupation takes place and 



charge localized in the system redistributes. Local charge 
is mostly accumulated on the upper level when applied 
bias value exceed £1 + Uu. Two-level system demon- 
strates such behaviour if the tunneling contact is sym- 
metrical (Figj3^i) or when system strongly coupled with 
tunneling contact lead k (sample) (FigJ3]3) . We have not 
found inverse occupation when two-level system mostly 
coupled with tunneling contact lead p (tip) (Figj3j;). In 
this case with the increasing of applied bias upper elec- 
tron state charge also increases but local charge continue 
being mostly accumulated on the lower electron state. 



We also analyzed tunneling current as a function of 
applied bias voltage for different level's positions (FigfSJ- 
FigHJi-f) . Tunneling current amplitudes presented in 
this work are normalized on 2r^ elsewhere. For all the 
values of the system parameters tunneling current de- 
pendence on applied bias has step structure. Height 
and length of the steps depend on the parameters of the 
tunneling contact (tunneling transfer rates and values of 
Coulomb energies). When both energy levels are above 
the Fermi level one can find six steps in tunneling current 
(FigJ2Ji-f). If both levels are situated below the Fermi 
level there are four steps in tunneling current (FigfJJi-f) 
and the upper electron level doesn't appear as a step in 
current-voltage characteristics but charge redistribution 
takes place due to Coulomb correlations. One can also 
reveal four steps in the case when only lower energy level 
is situated below the Fermi level (Fig|3Ji-f). 
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FIG. 4: Two-level system filling numbers a).-c). and tunneling current d).-f). as a function of applied bias voltage in the case 
when both energy levels are situated below the sample Fermi level. Parameters ei = —0.1, e-i = —0.3, U12 = 1.0, Un = 1.5, 
U22 = 1.6 are the same for all the figures. a),d).F fc = 0.01, F p = 0.01; b),e).F fc = 0.05, T p = 0.01; c),f).r fc = 0.01, T p = 0.03. 



IV. CONCLUSION 

Tunneling through the two-level system with strong 
coupling between localized electron states was analyzed 
by means of Heisenberg equations for localized states to- 
tal electron filling numbers taking into account high or- 
der correlations of local electron density. Various elec- 
tron levels location relative to the sample Fermi level in 
symmetric and asymmetric tunneling contact were inves- 
tigated. 



We revealed that charge redistribution between elec- 
tron states takes place in suggested model when both 
electron levels are situated above or below the sam- 
ple Fermi level. Charge redistribution is governed 
by Coulomb correlations. Moreover with variation of 
Coulomb interaction of localized electrons one can find 
the bias range of the two-level system inverse occupation 
when electron levels are localized on the opposite sites of 
the sample Fermi level. 

This work was partially supported by RFBR grants. 
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